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ABSTRACT 

Accreting white dwarfs (WDs) with non-degenerate companions are expected to 
emit in soft X-rays and the UV, if accreted H-rich material burns stably. They are 
an important component of the unresolved emission of elliptical galaxies, and their 
combined ionizing luminosity may significantly influence the optical line emission from 
warm ISM. In an earlier paper we modeled populations of accreting WDs, first gen¬ 
erating WD with main-sequence, Hertzsprung gap and red giant companions with 
the population synthesis code BSE, and then following their evolution with a grid of 
evolutionary tracks computed with MESA. Now we use these results to estimate the 
soft X-ray (0.3-0.7keV), H- and He II-ionizing luminosities of nuclear burning WDs 
and the number of super-soft X-ray sources for galaxies with different star formation 
histories. For the starburst case, these quantities peak at ~ 1 Gyr and decline by 
~ 1 — 3 orders of magnitude by the age of 10 Gyr. For stellar ages of ~ 10 Gyr, pre¬ 
dictions of our model are consistent with soft X-ray luminosities observed by Chandra 
in nearby elliptical galaxies and He II 4686A/H/3 line ratio measured in stacked SDSS 
spectra of retired galaxies, the latter characterising the strength and hardness of the 
UV radiation field. However, the soft X-ray luminosity and He II 4686A/H/3 ratio are 
significantly overpredicted for stellar ages of < 4 — 8 Gyr. We discuss various possibil¬ 
ities to resolve this discrepancy and tentatively conclude that it may be resolved by a 
modification of the typically used criteria of dynamically unstable mass loss for giant 
stars. 
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1 INTRODUCTION 

Accreting white dwarfs (WDs) are important for the 
study of binary evolution and accretion physics (see 
Postnov & Yungelson 2014; Maoz, Mannucci & Nelemans 
2014 for a recent review). Depending on the accretion rates, 
they are deemed to be components of cataclysmic binaries 
(see Warner 2003, for a comprehensive review) and super- 
soft X-ray sources (van den Heuvel et al. 1992). Under dif¬ 
fering initial conditions, they are the progenitors of com¬ 
ponents of future double white dwarfs, type la supernovae 
and binary millisecond pulsars (e.g. Whelan & Iben 1973; 
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Tutukov & Yungelson 1981; Webbink 1984; Iben & Tutukov 
1984; Taam & van den Heuvel 1986). Although they are 
widely studied, there remain many uncertainties concern¬ 
ing their formation and evolution. One avenue to constrain 
this problem is to invoke binary population synthesis. This 
allows one to compare the observed numbers and character¬ 
istics of these objects with that which are expected given 
our present understanding of binary evolution and accretion 
physics. 

Hydrogen-rich material accreted by a WD from its 
companion may burn at its surface either stably or unsta¬ 
bly. The result is strongly dependent on the mass trans¬ 
fer rate and mass of the WD (Paczynski & Zytkow 1978; 
Sion, Acierno & Tomczyk 1979). If the accretion rate cor- 
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responds to the stable nuclear burning regime, accret¬ 
ing white dwarfs (henceforth, SNBWD) have typical ef¬ 
fective temperatures of 10 5 — 10 6 K, emitting prominently 
in the soft X-ray and EUV. In this regime, not all ac¬ 
creting WDs can be observed as SSSs, as this depends 
on their intrinsic luminosity, effective temperature and the 
column density of the gas along the line of sight. In the 
present paper, we refer to SNBWDs with observed soft 
X-ray (0.30-0.70 keV) luminosity L x > 10 36 erg/s as 
SSSs. For accretion rates below the stability limit, hydro¬ 
gen burning occurs only in unstable flashes, observable as 
novae. There will be little X-ray emission in this regime, 
except in the post-novae SSS phase (Truran & Glasner 
(1995); Wolf et al. (2013), Soraisam et al. (in prep.)). 
If the accretion rate is larger than the maximum rate 
for stable nuclear-burning, the result is still unclear 
(Cassisi, Iben & Tornambe 1998; Hachisu, Kato & Nomoto 
1996). In the scenario of Cassisi, Iben & Tornambe (1998), 
the WD will become a red giant after accreting a small 
amount of mass (Paczyriski 1971). On the other hand, 
Hachisu, Kato & Nomoto (1996) proposed that in this 
regime the excess mass is lost as an optically thick wind. Fol¬ 
lowing Lepo & van Kerkwijk (2013), the WDs in this regime 
are dubbed rapidly accreting white dwarfs (RAWDs). In this 
scenario, an accreting WD’s photosphere will expand mod¬ 
estly, lowering their effective temperatures such that they 
radiate predominantly in the EUV. 

In addition to the emission from X-ray binaries domi¬ 
nating the X-ray radiation in most galaxies, extended emis¬ 
sion in soft X-ray band is widely observed in galaxies of 
different morphologies. However, the origin of the extended 
emission is still unclear. Bogdan & Gilfanov (2008) demon¬ 
strated that the emission from accreting WDs should be an 
important component of the unresolved emission. 

A population of accreting WDs with effective temper¬ 
atures T e ff > 10 S K should be capable of ionizing the in¬ 
terstellar medium. Rappaport et al. (1994) were the first to 
model ionization nebulae around SSSs, for the simple case of 
spherical symmetry and blackbody spectra. They found that 
these nebulae should have strong He II A4686, [O III] A5007, 
and [O I] A6300 emission lines. Woods & Gilfanov (2013) 
demonstrated that accreting WDs may provide a significant, 
or even dominant, contribution to the ionizing background, 
particularly in early-type galaxies with passively-evolving 
stellar populations (those with negligible ongoing star for¬ 
mation). Depending on their numbers, this could signifi¬ 
cantly impact the ionization state of the warm ISM in gas- 
rich early-type galaxies. Serra et al. (2012) found that many 
early type galaxies (at least 40% outside the Virgo cluster) 
host a detectable mass of HI, i.e. ~ (10 8 — 10 9 )Mq. Mor¬ 
phologically, this gas is found in regular HI discs, rings and 
irregularly distributed clouds. Most common are HI discs 
and rings, which may extend to tens of kpc, suggesting a 
covering factor of ~l/2 for the ionizing radiation from the 
stellar population. 

Although emission from populations of accreting WDs 
is important, it has been seldom studied. In Paper I 
(Chen et al. 2014), we conducted a population synthesis 
study of accreting white dwarfs. We found that it is impor¬ 
tant to have a realistic prescription for mass transfer and 
investigated the number of RAWDs, SNBWDs and the SNe 
la rate in the single degenerate (SD) scenario. In this paper, 



Figure 1 . Comparison of an accreting WD spectrum computed 
by means of a NLTE model (green solid line, Werner (1986, 1989); 
Rauch (2003)) with that found from the blackbody approximation 
(blue dashed line) for a WD with effective temperature T e fj = 
5.0 X 10 5 K, log(g/cm/s 2 ) = 8.0. The two vertical lines represent 
the H I and He II photoionizing limits. The edges around 25 A and 
15 A are the C VI (25.30A) and O VII (16.77A) absorption edges, 
respectively. The shaded area denotes the soft X-ray band (0.3- 
0.7keV) in our calculation. The NLTE spectrum was retrieved 
from TheoSSA (http://dc.g-vo.org/theossa). 


based on the results of our binary population synthesis of 
Paper I, we study the soft X-ray and UV emission of the 
population of accreting WDs and their contribution to the 
radiation of elliptical galaxies in the He II 4686A and H/l 
emission lines. These results are then compared with obser¬ 
vations. 

The paper is structured as follows. In section 2, we de¬ 
scribe how we calculate the emission spectra of accreting 
WDs. In section 3, we briefly outline the typical assump¬ 
tions underlying our population synthesis. In section 4, we 
present the soft X-ray (0.3-0.7keV) luminosity, H ionizing 
luminosity, and He II ionizing luminosity from the popu¬ 
lation of accreting WDs as a whole. We also compute the 
expected ratio He II4686A/H/3 as a function of age for a 
starburst stellar population and compare with observations 
of early type galaxies in the Sloan Digital Sky Survey (SDSS; 
Johansson et al. 2014). In section 5, we discuss the uncer¬ 
tainties in our standard model and present a potential solu¬ 
tion to resolve the discrepancy between the results found in 
our standard model and observations. In addition, we show 
the evolution of SSS number as a function of stellar age. 
Finally, we summarize our results and conclude in section 6. 


2 EMISSION SPECTRA OF ACCRETING 
WHITE DWARFS 

Detailed synthetic spectral models for accreting WDs can 
be produced using a non-local thermodynamic equilibrium 
(NLTE) stellar atmosphere model (e.g. Werner 1986, 1989; 
Rauch 2003), which is known to approximately reproduce 
SSSs spectra (e.g. Rauch et al. 2010). In Fig. 1, we compare 
the emission spectra of a SNBWD computed by means of 
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a NLTE model 1 with the blac.kbody approximation. Except 
for some absorption lines, the difference is not dramatic, 
particularly near the H and He II absorption edges. The 
difference in soft X-ray (0.3-0.7) keV luminosity between the 
blackbody approximation and NLTE models is smaller for 
accreting WDs with high temperatures and greater for low 
temperatures. With this discrepancy in mind, we adopt the 
blackbody approximation for SNBWD spectra in this paper. 

For a SNBWD, the nuclear burning luminosity is 

i'nuc = euXaMacc, (1) 

where eh = 6.4 x 10 18 erg/g is the nuclear energy release 
per unit mass of hydrogen, A'h is the mass fraction of hy¬ 
drogen and Mace is the accretion rate. We ignore the ac¬ 
cretion luminosity in the stable-burning regime, since it is 
much lower (less than 10% for Mwd < 1.2A7q and 20% 
for 1.20 Mq < Mwd < 1.4M©) compared with the nuclear 
burning luminosity. 

Detailed models of SNBWDs show that their photo- 
spheric radii are about 1 — 8 times as large as the radius 
of a cool WD (depending on the accretion rates, Iben 1982; 
Wolf et al. 2013). Wolf et al. (2013) simulated the evolution 
of accreting WDs with different accretion rates using the 
MESA code (Paxton et al. 2011, 2013). From their results, 
we know the dependence of the photospheric radius on the 
accretion rates in the stable-burning regime. For any given 
WD mass and accretion rate, we find the radius by linear in¬ 
terpolation of the results of Wolf et al. (2013) between WD 
mass and accretion rate. Then the effective temperature of 
SNBWDs is found from the Stefan-Boltzmann law. 

RAWDs with accretion rates larger than the maximum 
stable burning rate, expand as the WD loses mass and the 
effective temperature decreases. Hachisu, Kato & Nomoto 
(1999) investigated the optically-thick wind solution for dif¬ 
ferent white dwarfs with different accretion rates. They 
showed the dependence of photospheric radius and effec¬ 
tive temperature on the accretion rate (see their Figs. 3, 
4). By linear fitting of their results and linear interpo¬ 
lation among different WD masses, we can find the ra¬ 
dius and the effective temperature for any combination of 
white dwarf mass and accretion rate. It has been shown 
that RAWDs with typical mass accretion rates (less than 
4.0 x 10 _6 Mgyr _1 ) and effective temperatures (T « 1 
2 x 10 5 K) can fully ionize hydrogen and helium beyond the 
wind photosphere (Woods & Gilfanov 2013). For high accre¬ 
tion rates, the spectra of RAWDs are not well approximated 
by blackbody spectra because of absorption. However, high 
accretion rates can only be sustained for a short time and 
will contribute little to the total luminosity. Therefore, it 
is reasonable to use blackbody spectra as an approximation 
for RAWDs. 


3 BINARY POPULATION SYNTHESIS 

In Paper I, we have introduced our approach to studying 
populations of accreting WDs, and have shown the impor¬ 
tance of a careful treatment of the second mass transfer 
phase. Following the method outlined in Paper I, we carry 
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out two sets of calculations with our BSE+MESA model. 
Here we briefly summarize the procedure and main assump¬ 
tions of the calculation. 

First, we evolve a set of binaries with different initial 
parameters using the BSE code (Hurley, Pols & Tout 2000; 
Hurley, Tout & Pols 2002), and obtain the binary parame¬ 
ters at the beginning of mass transfer for a population of 
WD binaries with main-sequence, Hertzsprung gap and red 
giant companions. Second, we compute a WD binary evolu¬ 
tion library by means of the detailed stellar evolution code 
MESA (Paxton et al. 2011, 2013), which describes the evo¬ 
lution of WD binaries with different parameters at the be¬ 
ginning of mass transfer. For each WD binary with a given 
WD mass, donor mass and orbital period from the set of 
systems obtained in the first step of computations, we select 
the nearest track in the WD binary evolution library in or¬ 
der to follow its subsequent evolution. We take 1.2 Mq as 
the upper limit of the nascent CO WD, which was imple¬ 
mented in the BSE code, although modern estimates give 
« 1.07 Mq for single stars with Z = 0.02 (e.g. Doherty et al. 
2015). For the distribution of initial binary parameters, we 
use the Kroupa IMF for the primary mass (Kroupa 2001), a 
flat mass ratio distribution (Kraicheva et al. 1979) and a flat 
distribution in logarithmic space for binary separation (Abt 
1983). To describe common envelope evolution, we employ 
the Q-formalism (Webbink 1984; de Kool 1990). We use the 
fitting formula for the value of the binding energy parameter 
A from Loveridge, van der Sluys & Kalogera (2011). 

The efficiency a of the common envelope phase 
in ejecting the donor star’s envelope remains uncer¬ 
tain (see Ivanova et al. 2013, for a comprehensive re¬ 
view). The latest studies of post-common envelope bi¬ 
naries (e.g. Davis, Kolb & Knigge 2012) restrict a to a 
canonical value ^ 1 with a tendency to decrease with 
increasing secondary mass and WD mass. Some studies 
(e.g. Davis, Kolb & Willems (2010); Zorotovic et al. (2010); 
Ricker & Taam (2012)) show that, likely, a < 0.50. Here 
we adopt two values: a = 0.25 (model a025) and a = 0.50 
(model a050), performing calculations for each. In Paper I, 
we verified that the population of hydrogen-accreting WDs 
depends relatively weakly on the assumed value of a. 

Combining our binary population synthesis results for 
accreting white dwarfs with their predicted emission spectra, 
we study their X-ray and UV emission. As in Paper I, we in¬ 
vestigate two simple but representative cases. (I) Starburst: 
all stars are formed at t = 0, after which the star forma¬ 
tion rate is zero. (II) Constant star formation rate for 10 
Gyr. These two cases provide us with simple models which 
are representative of elliptical-like and spiral-like galaxies, 
respectively. 


4 RESULTS 

4.1 X-ray emission of accreting white dwarfs 

Figure 2 shows the evolution of the mass-normalized X- 
ray luminosity of the population of accreting WDs in the 
starburst and constant SFR cases. For our predicted ob¬ 
servable X-ray luminosity, we apply absorption with a col¬ 
umn density Ah = 3.0 x 10 20 cm -2 , which is a typical value 
for the extinction within the Galaxy at high Galactic lati¬ 
tudes (Dickey & Lockman 1990). In addition, for early-type 
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Table 1 . Comparison of the predicted X-ray luminosity in the soft X-ray band (0.3-0.7keV) in our standard model with observations. 
The observed X-ray luminosities for individual galaxies are taken from Bogdan & Gilfanov (2010); Zhang, Gilfanov & Bogdan (2012). 


Name 

age 

Lk 

A h 

L x observed 

L x predicted 

-^x,predict/ -^-'x, observed 

SSSs Num 


(Gyr) 

( l k,q ) 

(xlO^cm - ^) 

(erg/s) 

(erg/s) 



a) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

NGC3585 

3 

1.5 x 10 11 

5.6 

3.8 x 10 38 

(2.6 - 3.0) x 10 39 

6.8- 7.9 

330-570 

M32 

4-10 

8.5 x 10 s 

6.3 

1.5 x 10 36 

(1.2 - 1.7) x 10 36 

0.8- 1.1 

0 

NGC3377 

4 

2.0 x 10 10 

2.9 

4.7 x 10 37 

(2.0 - 2.2) x 10 38 

4.3-4.7 

30-60 

M31 bulge 

6-10 

3.7 x 10 10 

6.7 

6.3 x 10 37 

(3.9 - 6.1) x 10 37 

O 

1 

CO 

o 

10 

M105 

7-9 

4.1 x 10 10 

2.8 

8.3 x 10 37 

(1.7- 2.8) x 10 38 

2.0-3.4 

10-40 

NGC4278 

11 

5.5 x 10 10 

1.8 

1.5 x 10 38 

(0.88 - 1.98) x 10 38 

O 

05 

1 

CO 

10-40 


Notes. (l)-Galaxy name. (2)-Approximate stellar age of the galaxies. References- Terlevich & Forbes (2002); 
Coelho, Mendes de Oliveira & Cid Fernandes (2009); Olsen et al. (2006); Sanchez-Blazquez et al. (2006) (3) K-band luminos¬ 
ity of galaxies. (4)-Galactic absorption column density. (5)-Observed soft X-ray luminosity. (6)-Predicted soft X-ray luminosity. 
(7)-Predicted soft X-ray luminosity over observed soft X-ray luminosity. (8)-Predicted SSS number. 


galaxies we also show the observed X-ray luminosity per 
unit mass with column density 1.8 x 10“°cm~ 2 < Nh < 
6.7 x 10 2O cm _2 (the shaded region in the plot). This range of 
IVh corresponds to that observed for the six elliptical galax¬ 
ies in Table 1. In the calculation, we only considered Galactic 
absorption, since the intrinsic absorption of these elliptical 
galaxies is small. 

Our calculations predict that in the starburst case, the 
X-ray luminosity sharply decreases after 1 Gyr (Fig. 2). At 
early times (stellar age t < 1.0 Gyr), the X-ray luminosity in 
the starburst case is comparable with that in the constant 
SFR case for galaxies with the same mass. However, it is 
smaller for t > 1.0 Gyr. In order to compare our results with 
observations we use X-ray data from Bogdan & Gilfanov 
(2010) and Zhang, Gilfanov & Bogdan (2012) for nearby el¬ 
liptical galaxies. The data and our predictions are sum¬ 
marised in Table 1 where we list the observed X-ray lu¬ 
minosities along with the predicted X-ray luminosities and 
SSS numbers for six elliptical galaxies (Bogdan & Gilfanov 
2010; Zhang, Gilfanov & Bogdan 2012). Note that observed 
luminosities are total soft X-ray luminosities of unresolved 
emission and soft compact sources and therefore should be 
regarded as upper limits on the soft X-ray luminosity of nu¬ 
clear burning white dwarfs (see Bogdan & Gilfanov (2010) 
for details). In the calculation of the predicted X-ray lumi¬ 
nosities and SSS numbers, we have taken absorption into 
consideration for different galaxies with different column 
density values. 

One can see in Fig. 2 that, at the ages exceed¬ 
ing 6 Gyr, our predicted X-ray luminosities are compara¬ 
ble to that which are observed in nearby ellipticals, and 
is roughly two orders of magnitude below the soft X- 
ray flux expected if all SNe la are produced via the SD 
scenario (see Gilfanov & Bogdan (2010)). This is consis¬ 
tent with our prediction that the SD channel does not 
significantly contribute to the SN la rate at late de¬ 
lay times (see Fig. 8 in Paper I). However, for the two 
youngest galaxies in the Bogdan & Gilfanov (2010) sam¬ 
ple, NGC3585 and NGC3377, our standard calculations pre¬ 
dict up to an order of magnitude larger soft X-ray lu¬ 
minosities than observed. The strong dependence of the 


X-ray luminosity on the age of the stellar population 
predicted by our calculation seems to contradict obser¬ 
vations. In interpreting this result one should keep in 
mind that there is some controversy regarding the ages of 
the two youngest galaxies. Idiart, Silk & de Freitas Pacheco 
(2007) found that the age of NGC3377 is around 7.8 Gyr 
and Georgiev, Goudfrooij & Puzia (2012) listed literature 
values of the age of NGC3377, varying from 3.5 Gyr 
to 8.9 Gyr. In the case of NGC3585, Michard (2006) 
found a single stellar population age of 1.7 Gyr, while 
Hempel et al. (2007) suggested that its age should be 
larger than 3 Gyr. On the other hand, metallicity may 
be important to explain the discrepancy between model 
predictions and observations. O’Sullivan & Ponman (2004) 
and Georgiev, Goudfrooij & Puzia (2012) found that the 
metallicities of the two youngest galaxies, NGC3585 and 
NGC3377, are likely to be lower than solar metallic¬ 
ity. However, solar abundance is adopted in our calcula¬ 
tion. The metallicity plays a complicated role in the evo¬ 
lution of accreting WDs. The metallicity influences the 
initial-final mass relation for stars, the stable burning lim¬ 
its of accreting WDs, and the evolution of mass trans¬ 
fer rates (Umeda et al. 1999; Meng, Chen & Han 2008; 
de Mink, Pols & Yoon 2008; Kistler et al. 2013). Therefore, 
any strong conclusion based on the example of these two 
galaxies may be premature. On the other hand, a similar 
discrepancy was also found in comparing predicted popula¬ 
tions of low mass X-ray binaries with observations, based on 
a larger and different sample of galaxies and age determina¬ 
tions (Zhang, Gilfanov & Bogdan 2012; Fragos et al. 2013). 

4.2 UV emission of accreting white dwarfs 

Since nuclear-burning white dwarfs are characteristically 
high temperature sources, the emission from this popula¬ 
tion will ionize the ISM and will have an important influ¬ 
ence on the structure of the ionized gas and emission lines 
of a galaxy. In particular, compared with the single stel¬ 
lar population, accreting WDs are likely to emit predomi¬ 
nantly beyond the He II photoionizing limit (see Fig.l in 
Woods & Gilfanov 2013). So in this section, we will investi- 
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Figure 2. Mass-normalized X-ray luminosity (absorption ap¬ 
plied with Na = 3.0 X 10 20 cm -2 ) in soft (0.3-0.7keV) band 
for starburst case (upper panel) and constant SFR case (lower 
panel) as a function of stellar age. In the upper panel, the shaded 
area shows the X-ray luminosity with 1.8 X 10 20 cm —2 < Njj < 
6.7x 10 20 cm 2 . The blue solid and green dashed lines are for mod¬ 
els a025 and a050, respectively. The points with error bars are the 
observed X-ray luminosities for individual elliptical galaxies (see 
Table 1, Bogdan & Gilfanov 2010; Zhang, Gilfanov &; Bogdan 
2012). The cyan dashed line is the expected soft X-ray lumi¬ 
nosity if all SNe la are produced via the SD scenario assum¬ 
ing M = 10 _7 MQ/yr, and an initial WD mass of 1.2Mq and 
the delay time distribution given by Totani et al. (2008). The 
red dotted line is similar to the cyan dashed lines but assum¬ 
ing M = 3 X 10 _ 7 Mq /yr, an initial WD mass of l.OMg (cf. 
Gilfanov & Bogdan 2010). 


gate the H ionizing and He II ionizing luminosity of accret¬ 
ing, nuclear-burning WDs as a binary population. In addi¬ 
tion, we will compute the predicted ratio of two recombi¬ 
nation lines produced in the ISM, He II A4686A/H/3, and 
compare this with observations. 

In our calculation, we assumed the binary fraction to be 
50%. It has been suggested that the binary fraction may de¬ 
pend on the binary parameters (Kouwenhoven et al. 2009; 
Kraus & Hillenbrand 2009; Sana et al. 2012). We neglect 
this effect, and therefore may slightly underestimate the bi¬ 
nary fraction in our calculations. According to its standard 
definition, the binary fraction includes all binary stars, ir¬ 
respective of their separation and evolutionary state. Our 
population synthesis calculations keep track of only those 
binaries in which one of the components is a mass accreting 
white dwarf. Such mass-transferring binaries are in fact a 
very small minority in the population (in terms of the num¬ 
ber of systems and total mass). In accounting stars in wide 
binaries we assume that their emission is similar to the emis¬ 
sion of single stars of the same spectral type. To compute 
emission of single stars and stars in wide binaries we use 
the stellar population synthesis code of Bruzual & Chariot 
(2003), as described below. In computing the normalisation 
of the stellar emission we use the total mass of stars, ig¬ 
noring the small fraction (0.7% by mass ) of stars in mass 
transferring close binaries. Hereafter, we refer to the stars in 
wide binaries and single stars as single population (SP) and 
accreting, nuclear burning WDs as binary population (BP). 
The total emission of the entire population is computed as a 
sum of emission of the SP and BP. The SP-only model rep¬ 
resents the hypothetical case of absence of mass-accreting 
white dwarfs in the stellar population. 

In Fig. 3, we show the evolution of the H-ionizing and 
He-ionizing luminosities as a function of stellar age in the 
starburst case. For the binary population, we only show the 
results for model a025, since we do not find a significant dif¬ 
ference between models a025 and a050. To compute emis¬ 
sion from the SP, we use the Bruzual & Chariot (2003) stel¬ 
lar population synthesis model. In running their model, we 
assumed solar metallicity and Chabrier (2003) IMF, as this 
model is not available for Kroupa IMF. This introduces only 
small inconsistency with our BP calculations (conducted for 
the Kroupa IMF). The Chabrier IMF has the same expo¬ 
nent as the Kroupa IMF for stellar masses above 1.0 Mq 
and a more complicated shape at smaller masses. If the two 
IMFs are normalised to the same total mass, their high mass 
(> I.OMq) parts have the same shape and only small (~ 7%) 
difference in normalisation. This will result in a similarly 
small difference (~ 7%) in the predicted ionizing luminosity 
between the two IMFs, since stars with initial mass below 
1.0 Mq do not contribute to the ionising UV radiation. Such 
a small difference is insignificant given the large dynamical 
range of the ionising luminosities plotted in Fig. 3. In Fig. 
3, we do not show the case of a constant SFR, as in this 
case the H-ionizing and He-ionizing emission from accreting 
white dwarfs are significantly smaller than emission from 
the SP. This is due to the fact that in the constant SFR 
case there are many luminous young stars with high effec¬ 
tive temperatures, which dominate the ionizing emission. 

In the starburst case, the H-ionizing emission from the 
population of accreting WDs is comparable to that from the 
SP in the ~ 0.3 — 2.0 Gyr age range; it is unimportant for 
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Figure 3. Upper panel:H-ionizing (hv > 13.6 eV) luminosity per 
unit mass assuming ionization by single stars alone (SP only, blue 
dash-dotted line), ionizing radiation of a population of accreting 
white dwarfs for BP model a025 (a025 BP only, green dashed 
line), and their combined ionizing luminosity (green solid line) as 
a function of galaxy age in starburst case. Lower panel: similar to 
the upper panel, but for He-ionizing (hu > 54.4 eV) luminosity 
in starburst case. 

much earlier ages and for old galaxies. The He-ionizing emis¬ 
sion on the contrary is predicted to dominate the ionizing 
UV background at ages greater than ~ 0.3 Gyr. In the SP 
model, massive stars dominate the ionizing UV radiation at 
early times. As the galaxy age increases, the main sequence 
turn-off mass decreases and the effective temperatures of the 
hottest stars in the population drops. This leads to the de¬ 
crease of the ionizing luminosity seen in the Fig. 3 at t < 100 
Myr. Around t ~ 100 Myr, the first post-AGB stars appear, 
resulting in the sharp increase in the ionizing luminosity at 
t ~ 100 Myr, thereafter they dominate the UV emission. 
Note that in the SP, the He-ionizing luminosity provided by 
the massive stars drops very steeply and is outside the plot¬ 
ting range in the lower panel of Fig. 3. In the considered 
time range (t > 30 Myrs), the He-ionizing luminosity of the 
SP becomes significant only after the first post-AGB stars 
appear at t ~ 100 Myr. Note also, that the little “bump” 
around 1 Gyr, which can be seen on the curves in Fig. 3, 
is due to the fact that ~ 2 M© stars undergo a helium flash 
and their effective temperatures increase. 

Woods & Gilfanov (2013, 2014) proposed that the UV 
emission from accreting WDs is capable of ionizing the ISM. 
They investigated the contribution of accreting, nuclear- 
burning WDs in the SD scenario (as SNe la progenitors) 
and pAGB stars to the ionizing background and their influ¬ 


ence on the emission lines of warm ISM in early-type galax¬ 
ies. They predicted that, if the SD scenario for SNe la is 
the dominant channel, SNe la progenitors should strongly 
ionize He, producing strong He recombination lines in the 
extended emission-line regions of early-type galaxies. Such 
lines are predicted to be much weaker if only single pAGB 
stars are considered. Using a similar line of reasoning, one 
implication of the large population of accreting WDs found 
in our results across a large range of delay-times is that 
there should be strong He II recombination lines (e.g. He 
II 4686A) detectable in stellar populations with ages larger 
than 0.3 Gyr but less than ~ 5 Gyr. 

Among the recombination lines of He II, He II A4686 
is the strongest line in the optical band, and is not heavily 
absorbed by the ISM. In order to avoid the uncertainties 
in the covering factor, we directly compare the predicted 
ratio of He II 4686A/H/3 with observations. In order to 
predict the emission line fluxes, we make use of the pho¬ 
toionization code MAPPINGS III (e.g. Kewley et al. 2001; 
Groves, Dopita & Sutherland 2004). The procedure is sim¬ 
ilar to that in Woods & Gilfanov (2013). Here we briefly 
summarize the main points and assumptions in the calcula¬ 
tion. 

Serra et al. (2012) found that the common morphology 
of HI gas in elliptical galaxies is a regular HI disc or ring. 
The HI discs may extend to several kpc well beyond the 
nucleus, and may or may not be confined within the stel¬ 
lar body of the galaxy. Given that this gas will be ionized 
by the combined diffuse emission of many distant sources, 
we assume plane-parallel geometry in our calculation. We 
adopt a constant hydrogen density rh = 100 cm -3 , which 
agrees with the observed ratio of [SII] 6716A/6731A in 
passively-evolving galaxies (Yan & Blanton 2012). Regard¬ 
ing the metallicities of the gas, which can be estimated from 
the oxygen abundance, Athey & Bregman (2009) found an 
average value of Zoxygen — Z oxyg en, 0 . So we use solar abun¬ 
dances in our calculation. Although the hydrogen density 
may be varied as the galaxy evolves or the metallicity may 
be different, the line ratio calculated in this work is rela¬ 
tively insensitive to the density and the metallicity. Follow¬ 
ing Johansson et al. (2014), we assume that the ionization 
parameter log(U) ~ —3.5, where U = 47 r ^j p ^ ■ 

This is consistent with the [O III] A5007/H/3 ratio found in 
their stacks of early-type galaxies. 

In Fig. 4, we present the computed He II 4686A/H/3 
ratio powered by the emission from the SP together with 
the accreting WD population, and powered by the emission 
of the SP alone. At early times, the ratio begins to increase 
with the decrease in H-ionizing luminosity and increase of 
the He-ionizing luminosity. The ratio reaches its first maxi¬ 
mum when the H-ionizing luminosity reaches its minimum. 
Around 0.15 Gyr, the H-ionizing luminosity from the SP in¬ 
creases and exceeds the emission from the accreting WDs by 
an order of magnitude, leading to the beginning of a decrease 
in this ratio. After that, the H-ionizing luminosity does not 
change significantly. Thereafter, the evolution of this ratio 
simply reflects the evolution of the He ionizing luminosity. 
The luminosity ratio at 10 Gyr is consistent with observa¬ 
tions, but the computed ratio is larger than the observed one 
by a factor up to ~ 3 in the age range 1 — 8 Gyr. This comes 
about due to the predominance of accreting WDs in produc- 
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Figure 4. Comparison between the predicted values of He II 
A4686/H/3 in our starburst models and that which is observed in 
(stacks of) early-type galaxies, as a function of stellar age. The 
blue solid line shows the predicted values of He II A4686/H/3 for 
the combined populations (a025 BP+SP) and the green dashed 
line is for the single stellar population from (Bruzual & Chariot 
2003). The black dashed line shows the predicted values of 
He II A4686/H 0 for the model combining SNe la progenitors 
in SD-scenario and post-AGB stars (similar to the model in 
Woods &; Gilfanov (2013)). In the calculation of the emission of 
SNe la progenitors, we assumed that all the SNe la are produced 
via SD scenario, the initial WD mass I.IMq, WD effective tem¬ 
perature T e ff = 2 X 10 5 K and the delay time distribution given by 
Totani et al. (2008). The observed values (red squares) are data 
from Johansson et al. (2014). Note that for these points the verti¬ 
cal bars denote the error in the observed value, but the horizontal 
bars simply indicate the width of each age bin. 

ing our predicted ionizing background at these delay times. 
Although the total luminosities of accreting WDs and post- 
AGB stars in this age range are comparable, the lack of a 
sharp cutoff at 54.4eV in the spectra of accreting WDs leads 
to a dramatically greater production of He II-ionizing pho¬ 
tons. In particular, the He II-ionizing luminosity predicted 
in our a025 model exceeds that in the SP case by up to an 
order of magnitude at t~ lGyr (see Fig. 3). 

4.3 Emission from subsets of accreting WDs 

In Fig. 5, we show the soft X-ray emission from SNBWDs 
in model a025 and in an additional model in which we as¬ 
sume that upon exceeding the maximum accretion rate for 
stable burning, accreting WDs enter a CE phase instead 
of becoming RAWDs (model NORAWD) and WDs cease 
to contribute to X-ray and UV output of the population. 
Similar to Fig. 5, we show the H-ionizing luminosity and 
He-ionizing luminosity in Fig. 6. For the young stellar pop¬ 
ulations with age t < lGyr, the RAWDs dominate the H and 
He II ionizing emission. For old stellar populations (t > 1 
Gyr), the H and He ionizing luminosity is dominated by 
emission from the SNBWDs. This is due to the fact that, 
at early times, the binaries have massive donor stars and 
higher thermal timescale mass transfer rates. In contrast, 
at later times, the donors are less massive and the thermal 
mass transfer rate is lower, leading to more emission from 
the SNBWDs. In the model NORAWD, the soft X-ray lu¬ 
minosity, H-ionizing luminosity and He-ionizing luminosity 
are very different from that in model a025. In Paper I (see 


© 



Figure 5. Dependence of the soft X-ray luminosity (absorption 
applied with Ap = 3.0 X ! 0 20 cm 2 ) from different models of 
accreting WDs on the stellar ages in starburst case. The blue solid 
line is for the model a025. The green dotted line is for the model 
NORAWD assuming that accreting WD will enter CE instead of 
RAWD phases. 


Fig. 2), we found that there are two stable burning phases 
for some accreting WD binaries. In model NORAWD, re¬ 
placing the RAWD phase with a CE, only the first stably 
burning phase will contribute to the emission. So the H and 
He-ionizing luminosity in model NORAWD are comparable 
to the emission from SNBWDs at early and late times, when 
massive systems do not make a contribution. In Fig. 7, we 
show the computed Hell 4686A/H/3 ratio in the model NO¬ 
RAWD and in model a025. The ratio reduces significantly, 
but it is still larger than the observed ones by a factor of 2. 
This corresponds to a remaining discrepancy of up to factors 
of 2 to 5 between our prediction for the He II-ionizing lumi¬ 
nosity and that expected from the single stellar population 
alone (see Fig. 3). 

5 DISCUSSION 

5.1 Model uncertainties 

In the section above, we found that there is a strong de¬ 
pendence of soft X-ray luminosity To.3-o.7keV and the ratio 
He II 4686A/H/I on stellar age for non-star-forming early- 
type galaxies in our calculation, in contrast with observa¬ 
tions (Fig. 2 and Fig. 4). In producing a population synthe¬ 
sis model, there exist several uncertainties in both the initial 
conditions and the treatment of stellar evolution which may 
influence the results. In the following, we discuss a (partial) 
list of several aspects which may significantly impact our 
model, and their viability in providing a solution. Some of 
the model parameters not considered here were either tested 
and found to be of minimal importance (e.g. IMF, mass ra¬ 
tio), or are simply unlikely to be of significant impact (e.g. 
treatment of convective overshooting, as this only matters 
for massive stars). 

1) One important uncertainty is the retention effi¬ 
ciency of accreted matter, (see Bours, Toonen & Nelemans 






Hai-Liang Chen, T. E. Woods, L. R. Yungelson, M. Gilfanov and Zhanwen Han 



Figure 6. Upper panel: Evolution of H-ionizing luminosity from 
different subsets of accreting WDs as a function of stellar age 
in the starburst case. The blue solid line is for the SP model. 
The yellow solid line is for model a025. The green dash-dotted 
line and red dotted line are for the contribution of the SNBWDs 
and RAWDs, respectively. The black dashed line is for the model 
that assumes the accreting WDs will enter a CE instead of a 
RAWD phase. Lower panel: Similar to the upper panel, but for 
He-ionizing luminosity. 

(2013) 2 . In our calculation, we assumed that the He 
burning retention efficiency is 100%, which is clearly an 
overestimate and results in an excess of massive WDs 3 , 
which may be partially responsible for the high H and He 
ionizing emission of accreting WDs. On the other hand, 
Piersanti, Tornambe & Yungelson (2014) have shown that 
the conventional assumption that the retention efficiency of 
matter at the surface of a WD may be computed as a prod¬ 
uct of independently calculated retentions of H and He (i.e. 
rj = r]u x ?7He), results in an underestimate of r/. The rea¬ 
son is that H-burning heats the underlying He layer and the 

2 Note, in Yungelson (2010) He-retention efficiency is calcu¬ 
lated differently from that quoted by Bours, Toonen & Nelemans 
(2013). 

3 The WD mass may also be overestimated because of the initial- 
final mass function used in BSE code (see Claeys et al. 2014). 



Figure 7. The evolution of the HeII/H/3 ratio as a function of 
stellar age for different model assumptions in the starburst case. 
The blue solid line and red dotted line show the values for the 
combined population in model a025 and NORAWD, respectively. 
The green dashed line and square data are the same as Fig. 4. 



Figure 8. Similar to the upper panel of Fig. 2. The “new” model 
a025 is the result computed using the revised He burning retention 
efficiency (see text). 

outbursts of He-burning are milder than that in the case of 
accretion of pure He. 

In order to check if He retention efficiency has a 
strong influence on our result, we recomputed our re¬ 
sults with the He burning retention efficiency given in 
Hachisu, Kato & Nomoto (1999), and with the range in 
mass transfer rates corresponding to the stable burning 
regime as given in Wolf et al. (2013). In Fig. 8, we com¬ 
pare the soft X-ray luminosity from this new model with 
that from our standard model, and find the difference is not 
significant. In spite of this, we should emphasize that the 
retention efficiency suffers significant uncertainty, requiring 
further investigation. Therefore, without more precisely un¬ 
derstood retention efficiency, we can not completely exclude 
the influence of this parameter. 

2) In our calculations, solar abundance is adopted. As 
we discussed in Section 4.1, some galaxies in our X-ray com¬ 
parison are likely to have low metallicities, which plays a 
complicated and still uncertain role in the evolution of ac¬ 
creting WDs. On the one hand, at very low metallicities 
the stable-burning regime may be extended to much lower 
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mass transfer rates (Shen & Bildsten 2007). However, at the 
same time low-Z WD binaries are expected to have much 
higher accretion rates (Langer et al. 2000), likely pushing 
many more binaries into dynamically unstable mass trans¬ 
fer. Either of these mechanisms require very low metallicities 
(< O.OIZq), unlikely to be typical of many of the relatively 
nearby early-type galaxies used in our study. Therefore, we 
leave consideration of the evolution of WD binaries at low-Z 
to a future study. 

3) The X-ray comparison is based on a sample of six 
nearby galaxies (the only currently available set of measure¬ 
ments suitable for such a comparison). In addition, there 
is some controversy between different age estimates for the 
two youngest galaxies. However, while the incorrect age de¬ 
termination could in principle resolve the discrepancy in the 
X-ray band, it is unlikely to explain the He II 4686 A/ H/3 
line ratio, as demonstrated in Johansson et al. (2014). 

4) It may be that our modeling of the emission spectra 
of nuclear burning white dwarfs is incorrect. This must in¬ 
volve a failure of the ID theory, rather than any simple devi¬ 
ations from the black body spectral template used through¬ 
out the paper (see section 2). One obvious possibility is that 
the photospheric radii may be much larger than predicted by 
ID models of the nuclear burning on the WD surface. This 
would shift the peak of their emission to the UV or opti¬ 
cal bands where they would be remarkably unusual objects. 
However, Lepo & van Kerkwijk (2013) searched for such ob¬ 
jects in the Small Magellanic Cloud and found no plausible 
candidates. Therefore this possibility seems to be unlikely. 

5) One of the most important uncertainties in our 
calculations is the treatment of the CE phase. There are 
still many aspects of CE evolution which remain unclear, 
such as the available energy sources and the definition of 
the core-envelope boundary (see Ivanova et al. (2013) for a 
review, Hall & Tout (2014)). In addition, the criteria for the 
onset of a CE phase also suffer from great uncertainty. In 
the following subsection, we will discuss this possibility in 
depth and provide a possible solution to resolve the discrep¬ 
ancy between our computed results and observations. 


5.2 A potential solution 

In our calculations, for binaries harbouring giant donors, 
we adopted after Hjellming & Webbink (1987) and Webbink 
(1988) as the criterion for stable mass loss a critical value of 
mass ratio (hereafter, HW criterion): 


q c 


0.362 + 


1.0 

3.0(1 - rn c ) ’ 


( 2 ) 


where m c is the core mass fraction. For conservative mass 
transfer, this predicts dynamical instability for stars with a 
core mass fraction < 0.45, i.e. q c < 1.0. Chen & Han (2008) 
(see also Han et al. (2002)), using detailed binary evolution 
calculations, found that the critical mass ratio can be up 
to ~ 2.0, depending on the evolutionary phase of the donor 
star, mass and angular momentum loss from the system. 
Woods & Ivanova (2011) (see also Passy, Herwig & Paxton 
(2012)) found that this difference is due to the existence 
of a superadiabatic outer surface layer, which plays a crit¬ 
ical role during the mass transfer. Following on the work 
of Woods & Ivanova (2011), Pavlovskii & Ivanova (2015) 


found that the critical mass ratio varies from 1.5 to 2.2 for 
conservative mass transfer. Given these recent results, we 
turn now to considering alternatives to the HW prescription, 
in particular a more stringent condition for the onset of a 
CE, as a possible solution to the overproduction of accret¬ 
ing WDs found in our standard model above. We attempt 
to replace the HW criteria with q c = 1.5 (model a025qcl5), 
1.7 (model a025qcl7), and 1.9 (model a025qcl9), for bina¬ 
ries consisting of MS and giant stars. Although the mass 
ratio needed for the onset of a CE is understood to vary sig¬ 
nificantly for differing combinations of accretor mass, donor 
mass, and evolutionary state of the donor, these values are 
more consistent with the results of detailed numerical sim¬ 
ulations (e.g. Pavlovskii & Ivanova 2015), and much higher 
than what is typically found when implementing the HW 
criterion. In this section, using these new criteria, we calcu¬ 
lated the soft X-ray, H-ionizing and He-ionizing luminosity 
as a function of stellar age. 

In Fig. 9, we show the evolution of soft X-ray luminos¬ 
ity as a function of stellar age. Compared with our model 
using the HW criterion, the new criteria have larger critical 
mass ratios. With larger critical mass ratios, the predicted 
soft X-ray luminosity becomes lower. This is due to the fact 
that, with the requirement of a larger critical mass ratio, 
it will be relatively difficult for binaries to enter dynami¬ 
cally unstable mass transfer. Therefore, fewer binaries will 
enter a CE phase, and fewer accreting WDs will form. For 
old stellar populations, the predicted soft X-ray luminosi¬ 
ties in model a025qcl7 and model a025qcl9 are lower than 
the observed values. Given that the observational soft X-ray 
luminosity is the upper limit for the emission of accreting 
WDs, our results are consistent with observations. Fig. 10 
shows the evolution of the H-ionizing and He-ionizing lumi¬ 
nosity for different models, respectively. These results are 
similar to model a025, but the H-ionizing and He-ionizing 
luminosity is lower, for the same reason as above. In Fig. 11, 
we present the dependence of line ratio Hell A4686/H/3 on 
stellar age for the new models. For old stellar populations 
in model a025qcl7 and model a025qcl9, the line ratio is 
much smaller than that in model a025 and consistent with 
observations. 

The use of fixed unique critical mass ratios for the 
onset of rapid, unstable mass transfer is not realistic 
(Nelemans et al. 2000; Woods & Ivanova 2011). Numerical 
experiments above show that stability criteria of mass loss 
for giant stars play a crucial role in defining the number 
of SNBWD. However, stability also depends on mass of the 
donor, specific momentum of matter lost from the system 
etc. Fine-tuning of all these parameters may bring even bet¬ 
ter agreement of model and observations, but it is computer¬ 
time prohibitive. 


5.3 Number of SSSs 

As discussed in Paper I, not all SNBWDs can be detected 
as SSSs, since this is dependent on the luminosity and 
effective temperature of the source, and the column density 
of the gas along the line of sight. In this paper, we define 
SSSs as SNBWDs with soft X-ray (0.30-0.70keV) lumi¬ 
nosity L x > 10 36 erg/s. With this definition, we computed 
the evolution of the number of SSSs, which is shown in 
Fig. 12. In the starburst case, the SSS population peaks 
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Figure 9. Evolution of soft X-ray luminosity (absorption applied 
with Nh = 3.0 x 10 2O cm —2 ) as a function of stellar age in star- 
burst case. The red, green and yellow solid lines are for model 
a025qc!5, model a025qc!7 and model a025qc!9, respectively. 



Figure 10. Evolution of H-ionizing (upper panel) and He-ionizing 
(lower panel) luminosity as a function of stellar age in different 
models. 


around ~ 1 Gyr and then declines by ~ 2 orders of 
magnitude by the age of 10 Gyr. We note that the ionizing 
radiation from SNBWDs has a peak at the same age 
range (see the figures above). A similar peak was found 
in earlier studies of SNe la rates in the single degener¬ 
ate scenario (e.g. Canal, Ruiz-Lapuente & Burkert 1996; 
Han & Podsiadlowski 2004; Ruiter, Belczynski Sz Fryer 
2009; Yungelson 2010). The reason for this is that main- 



Figure 11. Evolution of line ratio Hell 4686/H/3 as a function 
of stellar age for combined population in model a025qcl5, model 
a025qcl7 and model a025qcl9. 




Figure 12. Evolution of the number of SSSs per unit stellar mass 
in the starburst case (upper panel) and the constant SFR case 
(lower panel). The blue dashed line is for model a025 and green 
dashed line for model a025qcl7 assuming Ah = 3.0 x 10 20 cm —2 . 
The red and cyan solid lines in the bottom panel show the number 
of SSSs with Ah = 3.0 x 10 21 cm —2 in model a025 and model 
a025qcl7. 


sequence stars with mass close to ~ 2 Mq have lifetimes 
~ 1 Gyr and radii and luminosities that enable a relatively 
long stage of thermal time-scale mass loss with the rate 
~ 10 _6 M©/yr, corresponding to the stable hydrogen burn- 
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ing regime (see Figs. 2 and 10 in Paper I) 4 . The number of 
SSSs for a 10 11 Mq galaxy at 10 Gyr is 10 — 20 in the star- 
burst case and 200 — 820 in the constant SFR case assuming 
Ah = 3 x 10 2O cm -2 . Given the high intrinsic absorption of 
spiral galaxies, we also made a set of calculations assuming 
Ah = 3 x 10 21 cm~ 2 . In this case, the SSS number is 40 — 130 
for a spiral-like galaxy. Note, however, that the observable 
number of SSSs depends on the absorption column density, 
which in turn depends on the relative distribution of gas 
and SSSs in any given galaxy. 

5.4 Gaseous nebulae around SSSs 

Given the hard emission of SSSs, it is suggested that they 
should be accompanied by ionized nebulae. Rappaport et al. 
(1994) modeled the ionization and temperature structure 
of such nebulae and predicted strong [OIII] A5007 and 
He II A4686 emission lines. Remillard, Rappaport & Macri 
(1995) searched for gaseous nebulae surrounding SSSs in the 
Large and Small Magellanic Clouds. They found only one 
SSS, Cal 83 (first identified as a high-excitation nebula by 
Pakull & Motch 1989), with a detected nebula, with null 
detections among the 9 other known Magellanic SSSs. This 
means that either the time-averaged luminosity and/or tem¬ 
perature of these sources must be much lower than presently 
observed, or that the local ISM density around SSSs are typ¬ 
ically too low. In fact, given that typical ISM densities are 
1-2 orders of magnitude below that observed for the Cal 
83 nebula, the latter appears more likely (see Woods & Gil- 
fanov (submitted), who found a 70% probability that < 1 
SSS in the LMC would have a detectable nebula). The situ¬ 
ation is complicated further in the accretion-wind scenario - 
for SSSs in star-forming galaxies, a wind-blown bubble will 
excavate a cavity in the ISM surrounding the source which 
is 10 - 40 pc in radius (Badenes et al. 2007). Consequently, 
the ionization parameter and morphology of these nebulae 
will be much different. 

Since the common morphology of HI gas in elliptical 
galaxies is a regular HI disk or ring, individual sources are 
not embedded in the ISM they ionize. With the above ex¬ 
planation and HI gas morphology in mind, we would ex¬ 
pect that the individual emission line luminosity in ellipti¬ 
cal galaxies should not be influenced by any complications 
introduced by the impact of the accreting WD on any sur¬ 
rounding ISM, as is the case in star-forming galaxies. 


6 SUMMARY AND CONCLUSIONS 

Based on the results of population synthesis of accreting 
WDs from Paper I and using simple assumptions regarding 
the emission of accreting WDs, we study the number of ob¬ 
served SSSs, the soft X-ray luminosity of accreting WDs, and 
their H and He II ionizing luminosity. The main conclusions 
are as follows. 

4 If mass loss is artificially fine-tuned, e.g., by the highly un¬ 
certain effect of mass-stripping effect of donors by an accretor 
wind (Hachisu, Kato & Nomoto 1999), the peak may be shifted to 
younger ages, e.g., Mennekens et al. (2010). The peak is smeared 
by the scatter in the actual parameters of WD+MS pairs - masses, 
mass ratios of components, and initial separations. 


(i) We compare predictions of our model with observations of 
nearby elliptical galaxies in optical and X-ray bands. To this 
end we use the He II 4686 A/ H/3 line ratio to characterise 
the hardness and strength of the ionizing UV continuum. In 
our standard model, using commonly (although not univer¬ 
sally) adopted assumptions, our predictions for the soft X- 
ray luminosity of old stellar populations are consistent with 
Chandra observations of several nearby elliptical galaxies. 
Likewise, for stellar ages of ~ 10 Gyr, the He II 4686 A/H (5 
line ratio from warm ISM predicted by our model is con¬ 
sistent with that measured in the stacked SDSS spectra of 
retired galaxies. However, for stellar ages of ~ 4 — 8 Gyrs our 
model significantly over-predicts the soft X-ray luminosity 
and the He II 4686 A/H/3 line ratio. 

(ii) We discuss various possibilities to explain this controversy 
and tentatively conclude that the most likely reason is that 
the classical Hjellming & Webbink (1987) criterion typically 
used for evaluating the onset of dynamically unstable mass 
transfer for giant stars predicts too low critical mass ratio. 
Replacement of this criterion by fixed critical mass ratios 
( q c = 1.7 or 1.9), brings the soft X-ray luminosity and line 
ratio He IIA4686/H/3 in the new models to consistency with 
the observations. However, we stress that the critical mass 
ratio needed for the onset of a common envelope is actually 
not a fixed value, and varies in an as yet unspecified way as a 
function of other stellar parameters. Further effort is needed 
in this regard, in order to better reconcile the predictions of 
binary population synthesis with observations. 

(iii) In the starburst case, the population of SNBWDs, their 
combined soft X-ray and UV output and the number of SSSs 
(SNBWDs with L x > 10 36 erg/s) peak at ~ 1 Gyr and then 
strongly decline with age, by ~ (1 — 3) orders of magnitude. 
Assuming Ah = 3-10 20 cm -2 , the number of SSSs present at 
10 Gyr in the starburst case is 10 — 20 for a galaxy with mass 
10 11 Mq. In the constant SFR case, assuming Ah = 3 • 10 21 
cm -2 , the SSSs number at 10 Gyr is around 40 — 130 for the 
same mass galaxy. 
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